Anomalous hydrodynamics and "normal" fluids in rapidly rotating BECs 
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In rapidly rotating bose systems we show that there is a region of anomalous hydrodynamics 
whilst the system is still condensed, which coincides with the mean field quantum Hall regime. An 
immediate consequence is the absence of a normal fluid in any conventional sense. However, even 
the superfluid hydrodynamics is not described by conventional Bernoulli and continuity equations. 
We show there are kinematic constraints which connect spatial variations of density and phase, that 
the positions of vortices are not the simplest description of the dynamics of such a fluid (despite 
their utility in describing the instantaneous state of the condensate) and that the most compact 
description allows solution of some illuminating examples of motion. We demonstrate, inter alia, 
a very simple relation between vortices and surface waves. We show the surface waves can form a 
"normal fluid" which absorbs energy and angular momentum from vortex motion in the trap. The 
time scale of this process is sensitive to the initial configuration of the vortices, which can lead to 
long-lived vortex patches - perhaps related to those observed at JILA. 

PACS numbers: 03.75.Kk, 73.43.-f 



C3 



-a 
o 

> 

&\ 
O 

in 
o 

-I— > 
c3 



i 

C 

O 

o 



X 



The area of rapidly rotating Bose Einstein condensates 
was one of the first to produce predicted phenomena 
quite distinct to the analogous condensed matter sys- 
tem, 4 He. Initially an instability was found 1] to a non- 
condensed, Laughlin, state for repulsive effective inter- 
actions between the atoms with sufficiently high angu- 
lar momentum, and the production^ of a fragmented 
condensate in the rotating attractive case. Both of these 
novel features occur in a regime where the atoms reside in 
restricted set of single particle states, the "Lowest Lan- 
dau level" (LLL), defined below. Subsequently it has 
been understood that when the atoms reside in the Low- 
est Landau level at intermediate amounts of angular mo- 
mentum mean field theory provides a good description 
(Mean field quantum Hall regime or MFQHR). In that 
regime the ground state is a vortex lattice^, Q] and it 
has been shown that other phases occur en route to the 
Laughlin state in the correlated domain[j| 0, Q as the 
amount of angular momentum in increased. The nature 
of these transitions in the thermodynamic limit and 
in finite systems remains a very active area of research, 
extending now into anisotropic traps|9lll0|. 

Following the pioneering studies UlI H^ ] of the vortex 
lattice, current experiments are in|l3j or near0,^| the 
mean field quantum Hall regime. The correlated regime 
is still to be investigated. The experiments find that 
the vortex lattice becomes very soft with an increasingly 
lengthy period for recovery after disruption. In this Let- 
ter we show the underlying hydrodynamics in the LLL 
is very unconventional and will be a contributory factor 
to understanding these experimental observations. There 
has been considerable theoretical discussion [1,^1 about 
the conventional or unconventional nature of excitations 
in the vortex lattice; this Letter provides a formulation 
of the hydrodynamics underlying such excitations. 



The defining characteristic of the wave functions, 
ip{x,y), belonging to the lowest Landau level is that 
ip{z,z*) = /(z)e~l z l / 2 , where z = x + iy and f(z) is 
analytic (or entire) in z. ip(z, z*) describes the system in 
the plane perpendicular to the axis of rotation. In gen- 
eral f(z, t) is time-dependent. The standard basis set, re- 
flecting the analyticity of f(z), is {(z m / y/irml) e~l 2 l Z 2 }, 
where m = 0, 1, • • •. 

To derive hydrodynamic equations from the Gross- 
Pitaevskii equation one conventionally represents the 
wave-function, tp = y/pe 1 ^, i.e. in terms of its modu- 
lus, y/p, and argument <f>. p is interpreted as the density 
and 4> is the order parameter, and in the nonrotating case 
is the velocity potential of the fluid. In the case of LLL 
hydrodynamics this motivates a representation of f(z,t) 
as f(z,t) = y^e"*, since the Gaussian factor is time- 
independent. The analyticity of / (i.e. df/dz* = 0) 
leads to "polar" analogues of the Cauchy-Riemann con- 
ditions for the derivatives of p and (f>: 
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The importance of these relations is that they con- 
strain kinematically the spatial variation of the density 
(via p = pe - ' 2 ' ) and phase unlike conventional super- 
fluid hydrodynamics. In addition, non-locality due to 
projection of the interaction term to the LLL affects both 
the continuity and "Bernoulli" equations, making them 
difficult to use. 

With these shortcomings in mind, we incorporate the 
analyticity of the wave- function from the start. There are 
clearly two complementary representations which 
accomplish this. Firstly the Taylor expansion of f(z, t) 
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may be written as: 

oo 

1>{z, z*,t) = f(z, t)e-^ 2 / 2 = a ™ (*) ^"-e-l-l 2 / 2 . 

m=0 

By the fundamental theorem of algebra, we may factorise 
the polynomial in terms of its roots Ca(t) %$)■ 

oo 

where S is the normalisation integral. The points Ca{t) = 
X a {t) + iY a (t) are the vortex positions, or coordinates, 
at time t. 

We immediately see that there is an important dis- 
tinction between a condensate wave-function residing in 
the LLL and a general condensate wave-function. In the 
former case the wave-function is completely specified by 
the position of the vortices. In particular, the density is 
specified by the vortex positions - there are no indepen- 
dent sound waves (long wavelength Bogoliubov quasipar- 
ticles), unlike the general case. There is thus no possi- 
bility of the existence of a normal fluid component in a 
conventional sense. 

Since the vortices uniquely determine the wave- 
function in the LLL, it is natural to re-express the Hamil- 
tonian in terms of the vortex coordinates, C, a (t). We use 
the standard model Hamiltonian for N particles in a ro- 
tating reference frame: 

N N N N 

^ = -iE v «+5E^+^ E *(r n -rv)-^E L « 

n— 1 n— 1 n^n' — l n—1 

Units of length, £q, and energy, hu>o, are those provided 
by the harmonic trap; angular momenta, L z n are scaled 
by h. There are two remaining dimcnsionless parameters. 
Firstly, oj, is the angular velocity of the rotating frame 
divided by the natural frequency of the harmonic trap. 
Secondly the coupling constant, rj = Annad"^. 

Then the Hamiltonian takes the form (assuming M 
vortices): V, inter-particle potential, U energies and the 
rotating frame term, 1Z, are defined as (using ' = d/dt): 

M 

S = nNY / \P*M- m (0\ 2 ™<- 

m 

1 M 

IC = V = -ttJV^ \P* M - m (0\ 2 (™ + 1)! 

~ m 

M 

n = nNcjJ2\P*M- m (0\ 2 mml 

rn 

U = \N{N - 1)S- 2 E™,„, P , 9= P*M-rn(OP*M- n (0 

P M - P (C)PM- q (C)(p + q)V- (p+q) 5 m+n , P+q 



where the quantity P n (Q = £ il<i2 <... <in CiiC< 2 ' is 
the nth symmetric polynomial in the variables Co- 
lt is instructive to contrast this expression with the 
corresponding expression for vortices in an incompress- 
ible medium in a container of size R: 

M 

H = -\T 2 n Q J2 ln \Ca-Cp\-^MP 2 -\U 2 )+ images 

where u> is the unsealed angular velocity of the rotating 
frame, L is the circulation of the vortices (assumed to be 
identical) and uq is the constant density. 

There are two striking differences. Firstly the in- 
teraction between the vortices is analytic and is not pair- 
wise, with multi (up to M-) vortex contributions. Sec- 
ondly the expression for the angular momentum of the 
fluid in terms of the vortices (entering the term involv- 
ing u>) is not a sum of single vortex contributions but is 
collective in nature in the LLL case. 

To determine the equations of motion of the conden- 
sate constrained to reside in the LLL, we employ a vari- 
ational formulation^^ El where the LLL wavefunction 
is a trial function leading to the Lagrangian, C = T — 7i, 
where (assuming ip(z,z*,t) is normalised) 

M 

T = nNY / rn\^(P^_ m (OPM- m (0-PM-m(OPM- m (0) 

m=0 

If the corresponding action is varied with respect to the 
( a the resulting equations of motion for the vortex co- 
ordinates are rather complicated. However the form of 
C shows that the vortex coordinates always enter in the 
form of the symmetric polynomials, P n ({Ca})- This sug- 
gests that their numerical values should be used as the 
primary variables. Since, up to a sign, the P n ({Ca}) are 
the coefficients in the polynomial whose roots are the ( a , 
there is a one-to-one relation between the vortex coor- 
dinates and the symmetric polynomial values. The dy- 
namics are most simply expressed in the variables p m (t) 
and m (t) implicitly defined by: 

(-irP U -m(0 =a m = x [^e w ™ 

V 7rm! 

In particular normalisation can be easily enforced via the 

Pm- 

In terms of p m and 9 m we find ('indicates a time deriva- 
tive and A = X(N - 1)/4tt): 

{M 
S' 1 E Pml-Om ~ (1 + m[l - to])} 
m=0 

M 

- A E ^ PmPnP pPq ^ e ^-^ ] 
m,n,p,q—l 

I (m + n)! / (p + g)l s \ 
V 2™+«m!n!Y 2P+ip!q\ p+q > m+n j 
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where S = E m =o Pm- 

Differences between incompressible]^ or Thomas- 
Fermi j2^| and LLL vortex dynamics are already appar- 
ent in the case of m vortices symmetrically arranged and 
equidistant with respect to the centre of the trap. Up to a 
phase, the vortices are in the directions of the mth roots 
of unity in the £-plane, and hence the wave-function is 
of the form <j)(z) ~ (ao(t) + a m (t)z m ). We will initially 
restrict ourselves to a reduced Lagrangian only involving 
the p m , Oq and 9 m variables (as po can be eliminated by 
the normalisation constraint 1 = po + p m ): 



N 



-Oo - 1 - A 

Pm ([-d m - Oo] - m[l -co}- A[2-(™- 2 ) - 2] 
+pl~\(2-( m -V _i_1^2- 2 " ; 



\\2- 



(ml) 



8q is not determined and its time derivative may be in- 
terpreted as the zero of energy in the original quantum 
problem. 

The significant Euler-Lagrange equation is that associ- 
ated with p m (the equation for 6 m expresses the unitarity 
of the original time evolution, in this case p m = 0) , yield- 
ing the excitation frequency uj m : 

uj m = —6 m + 9q 

= m(l-^)-2A(l-2- (m ' 1) ) 



-p m 2X 



(ml) 2 



This expression describes, in a unified manner, several 
physical phenomena: surface wave excitations, instabili- 
ties and anomalous vortex dynamics. 

We may re-express the equation for u> m in terms of 
the radial distance of the vortices from the centre, £, as 
p m = 1/(1 + £ 2m /m!). This implies a natural scale for 



the radial vortex positions of (,„ 



(ml) 



|U/(2m) 



m, 

which accords with the extent of the condensate with m 
quanta of angular momentum. For £ < Q rni the vortices 
are "inside" the condensate, and for ( > Q m the vortices 
are in the evanescent tail. 

Thus in the limit that the vortices are very far from 
the centre of the trap (i.e. p m — > 0) we recover the linear 
surface excitation frequencies derived by Kavoulakis et al 
|l9| . We are thus immediately provided with a represen- 
tation of surface waves in the LLL of angular momentum 
m in terms of m vortices. A relationship between vor- 
tices and surface waves in the Thomas-Fermi limit has 
been discussed by Anglin [2(| and Tsubota et al|2l|. 

To obtain the equilibrium radial distance of the vor- 
tices, Ceq(w)j w e set uj m = so that the vortices are 
stationary in the rotating frame. We find Ceq(^) = °o 



for lo less than a critical rotation frequency, to, 



(m) 



1 - 2X 1 - 



tor to > ujc 1 we see: 



Cc 2 q m M 



2A 



/ (m) \ 

m(u> — ujc ) 



(2m)! _ 2 
(ml) 
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These critical frequencies do not have a simple physical 

interpretation as shape instabilities of the condensate - 

(2) 

as, for example, coincides with the frequency of a 
first order transition of a single vortex to the centre of 
the trap. This is unlike the Thomas- Fermi limit . 

One might expect the dynamics of the vortices when 
they are "inside" the trap, £ < £ m , to be reminiscent 
to incompressible vortex dynamics as the phases of the 
wave-functions are identical. However the LLL dynamics 
is completely different. For example as the separation of 
the vortices vanishes, the frequency Lu m {C — 0) obeys: 



^m(C = 0) 



m(l-u) + 2X2 



-Ira 



(2m)! 
(m!) 2 



m(l — lo) + 



2A 



It does not diverge, in stark contrast with the incompress- 
ible case. There the frequency, O, of rotation of such a 
polygonal array of m vortices, each with circulation T, a 
distance £ from the centre is|22j (in the laboratory frame) 
n = T(m- 1)/(4ttC 2 ). 

We will now show that energy and angular momentum 
can be deposited in surface waves (vortices at large dis- 
tances) by vortices near the centre of the trap. In other 
words the surface waves can in part play the role of a 
"normal fluid" . Consider n vortices which are all initially 
very near the centre of the trap. A linearised approxi- 
mation (where p n ~ 1) allows coupling to all the angular 
momentum channels which conserve angular momentum 
(i.e. m < 2n). This implies we will have n additional vor- 
tices, which will turn out to be at large distances from the 
centre of the trap. Because of angular momentum con- 
servation p n — m — Pn+m is a constant, which we choose to 
be zero. Then the equations for R m = \(p n -m + Pn+m) 



and ov 



(0n- 



1 \ u ri+ra "nf c 

—A m R m sincr„, 
Xm - A m cos a Tl 



Where 
A m = 4A- 



f2n)! 



2 2 "n! (n + m)l(n — mjl 



X:, 



= 4A 



2 2n n\ 



(2n) 



(2n 



(2n- 



l (n — m)! 



l (n + m)l 



These equations are readily integrated and one finds ex- 



ponentially unstable growth for n + Cy 



> 



m > n + l 



and n — l>m^n — c-y/n, in the limit where n 3> 1, 
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where c = 2(— ln[V2 — l]) 1 ^ 2 - The extent of this un- 
stable region is determined by the spatial overlap of the 
states with the state n. Oscillatory behaviour occurs for 
other values of m. The time constant, r m , for the growth 
is t,^ 1 = | \J A^j — Xm ■ For the most divergent cases, 

m* ~ n±2(ln[2/(V5-l)]) 1 / 2 ^, we find r^ 1 « A/^ 
If we assume that initial values of p m , p m (0), are of 
the order of e <C 1, and then the unstable modes grow 
by time t to p m * ~ 6 where e <C 5 <C 1, then we can 
deduce the corresponding vortex positions from the crude 
representation of the polynomial: 

z" + 5{z n -^ + ■ ■ ■ + z n ~ l + z n+1 + • • • + z n+ ^) 

+ e(l + z + ■ ■ + z n -^- 1 + z »+v^+i + . . = o 

Noting that if 2 is a solution so is 1/z, let us look for solu- 
tions for large \z\. Keeping the largest terms multiplying 
each of the two small parameters, we obtain 

1 + x + -^x^ ~ 

where z^™ = x/5. Balancing the middle term with either 
the final term or with unity, we find y/n approximate 
solutions z'v™ = —1/5 and n — y/n approximate solutions 

^Vn(Vn-l) = _£/ £ _ ( We assume e < §y/n.) 

Geometrically, the roots form four concentric rings, 
two inside the trap and two outside. The closest and 
the most remote rings do not change their moduli signif- 
icantly, but the rings with moduli l^ 1 /^/™ change with 
time. The inner ring of y/n vortices, moving outwards 
and the outer ring (with the same number) moving in- 
wards. Each member of the inner ring has a correspond- 
ing member of the outer ring with the same argument. 
This demonstrates the exchange of angular momentum 
between the inner and outer rings. 

This description assumes that one draws the initial 
conditions from the same probability distribution for 
each coefficient in the polynomial - which leads to rings of 
roots [25I Hi| . However for initial vortex positions form- 
ing a cloud distributed normally around the origin in the 
complex plane, then the coefficients of the polynomial 
are very strongly peaked [2f| at m ~ n/2. Such a config- 
uration will evolve much more slowly as it has negligible 
initial weight in the unstable modes, n — y/n < m < n. A 
simple estimate using the expressions above yields time 
scales for evolution i r ; ng ~ y/nlnn and t c i ouc i ~ r?l 2 . It 
is tempting to compare this distinction with the JILA 
experiments |27j where long-lived concentrations of vor- 
tices were observed. Although the experiments were not 
conducted under LLL conditions on average, in the low- 
density environment of the vicinity of the vortices the 
dynamics might be of a LLL nature. 
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